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We bound the coupling of pseudo-scalar particles to TtGAGqcd using (the lack of) monojet plus 
missing Et events at the Tevatron, and estimate the bounds obtainable from LHC. In addition, we 
revisit the bounds on the coupling to FA -Fem from e + e~ collider events with single photon and 
missing Et final states. This is especially interesting in light of the recent experimental results from 
the PVLAS collaboration, which we believe can be tested by data which will be available in the 
near future. 
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INTRODUCTION 

In this note we use collider data to bound the cou- 
plings of pseudoscalar fields to gluons and photons. The 



processes we consider are e + e 
liders, and pp or pp — > single jet 



7+ $i 



col- 



Jfjrp for hadron collid- 



ers. The signature is similar to that for gravitons flow- 
ing into extra dimensions, although they can be distin- 
guished. For instance, one characteristic of the pseudo- 
scalar amplitudes is that the cross section is indepen- 
dent of center of mass energy at high energies. This fol- 
lows from dimensional analysis: the interaction g<pF f\F, 
where <f> is the pseudoscalar field, has a coupling constant 
g with dimensions of inverse mass. Therefore 2 — > 2 pro- 
cesses involving the production of a single (f> particle will 
have a cross-section da/dVL = g 2 f(s/t), where f(s/t) is a 
function which depends only on the angle. 

Specifically, the couplings we will bound are the fol- 
lowing [lfij : 



for gluons, and 



(0.1) 



(0.2) 



for photons. We will ignore any other possible couplings. 

There are a number of motivations, both theoretical 
and experimental, to be interested in such interactions. 
Recently, the PVLAS collaboration has reported data 
consistent with the existence of a light, neutral, pseu- 
doscalar particle coupled to the photon as above |lj . The 
favored mass range is 0.7 meV < < 2.0 meV, and 
the coupling 1.7 lO" 6 GoV" 1 < g < 1.0 10~ 5 GcV" 1 . 
These values appear to be in contradiction with astro- 
physical bounds coming from stellar dynamics and the 
lack of detection of such particles produced in the sun 
(see section 4. of for a recent review); however in 
some models these constraints might be relaxed In 
light of this, it is useful to have a direct bound on g com- 
ing from particle physics data. As we will see, an analysis 
of currently available data would come to the edge of this 
region of parameter space, and in the near future data 
will be available which will constrain it. 



One bound on the coupling g comes from decays of 
the T meson to photon plus missing energy. The exper- 
iment CLEO Q bounds the branching ratio T(IS') — > 
<j>j < 1.3 10~ 5 , which translates (assuming the </> cou- 
ples only to photons and not directly to quarks) into g < 
1.0 10~ 3 GcV -1 . A better bound g < 5.5 10" 4 GeV" 1 
was obtained using e + e~ collider data from ASP Q Q. 
As we will see, current collider data could improve this 
by nearly two orders of magnitude. 

One particularly interesting pseudoscalar particle is 
the QCD axion, proposed as a solution to the strong 
CP problem and which can also serve as a dark 
matter candidate. For our purposes, there are two pri- 
mary differences between axions and more general pseu- 
doscalar particles with couplings (|0.1|) and (|0.2|) . The 
first is the relation between the axion coupling and mass: 
m a ~ AQ CD /f a . For bounds obtainable from colliders, 
the implication is that the mass of the axion can be effec- 
tively set to zero in the parameter range we are consid- 
ering. The second is the relation between f a and ga- 
ga = (aC) / (2n fa) , where C is a model-dependent con- 
stant which is typically of order 1, but in special models 
can be significantly smaller. 



SUMMARY OF RESULTS 

In this section we will summarize our estimates of the 
bounds on the pseudoscalar couplings which could be 
obtained with an analysis of collider data. The details 
of our analysis can be found in the sections that follow. 
We will give all bounds at the 95% confidence level, and 
for experiments where the data have not yet been ana- 
lyzed we will quote bounds under the assumption that the 
results agree with the central value for the standard model 
background. 

Since we are looking for a missing energy signature, we 
will require that the pseudoscalar does not decay inside 
the detector. Therefore the efficiency should include a 
factor of exp(— LT/j), where T is the inverse lifetime, 7 
is E^/rrifp, and L ~ 10 m is the size of the detector. The 



2 



partial width for decay to photons is 



■ <A^7 7 



2 3 

9 m% 
64tt ' 



(0.3) 



For the case of non-zero / and zero tree-level g, if the 
axion is too light to decay hadronically, the dominant 
decay mode will be to photons via a two-loop diagram 
involving gluons and charged hadronic matter. With- 
out assuming a specific model we can not determine pre- 
cisely what the width will be; however generically this 
rate is small and the will not decay in the detector. On 
the other hand, if is large compared to the mass of 
the pion the axion can decay hadronically and it is diffi- 
cult to make a model-independent estimate of the width. 
Therefore we will take the bound on the axion mass to 
be nicf, ~ m„o. 

In a particular model with calculable branching ratios, 
the bound at larger masses could be improved by look- 
ing for displaced vertices in the detector. However, in 
what follows we will simply quote mass bounds using the 
requirement that the <fr travel for a distance greater than 
ten meters, and assuming that the coupling saturates the 
stated bound. 

Using the data from the D0 experiment at Fcrmilab: 



/ > 35 GeV. 



(0.4) 



A bound on the photon coupling g was obtained in 
for the process e + e" — > 7+ missing energy: g < 
5.5 10 _4 GcV _1 for a pseudoscalar mass < 25 MevQ~. 
This bound comes from the ASP detector at the Stanford 
Linear Accelerator Center storage ring PEP. 

We have analyzed several past and future experiments 
to determine the bound which could be obtained assum- 
ing the data are consistent with the standard model back- 
grounds. Since the amplitudes are independent of energy, 
the bounds can be improved mainly by increasing the 
total luminosity. An analysis of the combined data from 
LEP at ALEPH, OPAL, L3 and DELPHI would yield the 
more restrictive bound: 



g < 1.5 10~ 4 GeV -1 



(0.5) 



for m < 65 MeV. 

As for future and ongoing experiments, we expect the 
Large Hadron Collider at CERN to improve the bound 
on /: 



/lhc > 1300 GeV. 



(0.6) 



For the PEP-II e + e collider, current integrated lumi- 
nosity gives 



5PEPH < 8.9 10~ 6 GeV -1 
for m A < 0.12 GeV. 



(0.7) 
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FIG. 1: 95 % confidence bounds on pseudoscalar-gluon cou- 
plings versus pseudoscalar mass obtainable for hadron collid- 
ers. The testable regions for each experiment are below the 
line. 



A similar bound can be obtained from KEKB e + e~ col- 
lider. With the current integrated luminosity the bound 
would be 



<?kekb < 8.2 10~ 6 GeV" 



(0.8) 



for < 0.13 GeV. The expected total luminosity for 
KEKB is at least twice the current total, which would 
improve the bound to g < 5.9 10~ 6 GeV _1 . 

Finally, for the Super KEKB upgrade to KEKB is 
expected to produce 10 7 pb^ 1 per year. After two years, 
this would improve the bound to 



5kekb<1.9 HT 6 GeV" 



(0.9) 



which would rule out most of the parameter space favored 
by PVLAS. 

In figures n an d |21 we summarize the exclusion regions 
for these colliders. 



BOUNDS FROM HADRON COLLIDERS 

In this section we analyze the bounds on the pseu- 
doscalar gluon coupling that can be obtained from pp 
annihilation to one jet plus $ T . The relevant subpro- 
cesses are represented in figure |3J 

The data obtained by the D0 detector at the Teva- 
tron between 1994 and 1996 corresponds to an integrated 
luminosity of 78.8 ±3.9 pb _1 with a center of mass energy 
= 1.8TeV. The backgrounds for this process were 
analyzed in [8j, in which the authors obtained bounds 
on a similar process in which the missing energy was 
lost into an extra dimension via Kaluza-Klein graviton 
emission. For our analysis we will use the same cuts 
as 0: $ T > 150 GeV, E T > 150 GeV, and pseudo- 
rapidity < 1.0. The dominant background is the 
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FIG. 4: Feynman diagram of the process e e — * 'yep. 

section for this background is a u p ~ 4 10~ 3 pb, which 
corresponds to 400 events. If we estimate 5% error in the 
background and require 95 % confidence (N^ ~ 40), we 
need er^ ~ 4 10~ 4 pb. 

Including the cut, this corresponds to a bound on / 

/ > 1300 GeV. (0.11) 



FIG. 2: Bounds on pseudoscalar-photon coupling versus pseu- 
doscalar mass for e + e" colliders, the preferred region for the 
PVLAS results, and the band covered by typical QCD axion 
models. The area above the lines could be tested with an 
analysis of collider data. 






FIG. 3: Feynman diagrams for the process pp 



' + jet. 



production of a Z which decays to vv, where the neutri- 
nos carry the missing energy. The results of a numer- 
ical simulation using the software package CompHEP 
are that the cross-section for monojet plus axion is 
<Jpp->^t = 293 pb (GeV//) 2 . 

Using the error estimates of Q and requiring 95% 
confidence, the bound obtained from this cross-section 
(assuming a detection efficiency of 1 given these cuts) is: 



/ > 35 GeV. 



(0.10) 



A similar analysis could be performed using the data 
from CDF • The luminosity considered in the analysis 
of 0| is 84p6 _1 . We expect the result to be similar. 

As for the Large Hadron Collider at CERN, we will 
take an integrated luminosity of 10 5 pb -1 , and select 
monojet events with £r,jet > 1 TeV. As we will see, given 
our estimates the signal will be limited by systematic 
error in the background rather than by statistics. Again, 
the primary standard model background will be pp — > vv 
+ jet. At center of mass energy y/s = 14 TeV, the cross- 



BOUNDS FROM e+e COLLIDERS 



Consider the process e + e~ — > j<j) (see figure 0J. The 
differential cross section is (in the limit yfs » m e , but 
including the mass of the <f> particle): 



d^J cm 



ag 



(s - m 2 ) 



1287T 



(0.12) 



and the total cross section is ot 



ag 2 («-"4) 2 
24 T 2 ' 

One relevant fact is that the signal we are looking for 
has a two-body final state, with both particles effectively 
massless. As a result, the photon always carries precisely 
half the center of mass energy. On the other hand the 
standard model backgrounds have at least a three body 
final state, with zero phase space for the photon to carry 
this energy. Furthermore, if the masses of some final state 
particles are significant, there is a further suppresion due 
to conservation of energy and momentum. If the energy 
resolution of the detector were infinite, we could therefore 
cut out all the standard model background by requiring 
the photon energy to be very close to the beam energy. 
In a real experiment, uncertainties in the photon energy 
will be relevant, and we take that into account in what 
follows. 

As we mentioned above, this process was considered in 
0, using data from the ASP detector at SLAG Improved 
bounds can be obtained using LEP data from ALEPH, 
OPAUL3 and DELPHI. For instance, the ALEPH detec- 
tor [ll| took data in a range of center of mass energies, 
between 188 and 209 GeV, with integrated luminosity 
628 pb" 1 . The collaboration performed an analysis of 
single photon events, and found that the data are con- 
sistent with the predicted standard model background, 
which is dominated by "/vv in the final state. However, 
for reasons discussed above a better bound on g could be 
obtained with more stringent cuts on the energy of the 
photon. The energy resolution of the ALEPH calorime- 
ters is of order 3 GeV for a 100 GeV photon. For example, 
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for a beam energy of 94 GeV, if we require E 1 > 91 GeV 
and |cos(0)| < .96, o e ^ lvJ> = 2.7 10~ 3 pb, compared 
to cr e<! ^ 7 = (g/10- 5 ) 2 1.2 10~ 5 pb for the signal. The 
expected number of back grou nd events is then 2, and 
requiring 95% confidence [15| (4.72 events) would give 
a bound (assuming the results agree with the standard 
model) 

g < 2.6 10~ 4 GeV -1 , (0.13) 

valid for < 48 MeV. 

Similar results can be obtained from the other LEP 
experiments: L3 0, with 189 < y/s < 209 and lumi- 
nosity is 619 pb" 1 ; OPAL 0, with y/s = 183 GeV, 
and luminosity is 54 pb -1 and DELPHI 01, with y/s = 
130 — 209 GeV and luminosity is 667 pb -1 . Combin- 
ing these gives roughly a factor of four in luminosity, but 
with these cuts we must require 6.8 events for 95% confi- 
dence, and therefore the bound is g < 1.5 10~ 4 GeV -1 , 
for m < 65 MeV. 

Better bounds can be obtained by higher luminosity 
colliders: PEP-II and KEKB. As explained above, requir- 
ing strict cuts on the energy of the photon (in the center 
of mass frame) greatly reduces the backgrounds due to 
many particles in the final state and particularly mas- 
sive particles (BB processes, K^KlJ processes where 
the Kl's are interpreted as missing energy) and we do 
not expect them to be significant. The energy resolu- 
tion of the BABAR and BELLE detectors for a 5 GeV 
photon (in the center of mass frame) is around 1.5-2 %. 
Therefore in the following analysis we will require that 
E 1 > 0.985E BEAM . 

The center of mass energy of both colliders is sim- 
ilar, as both are tuned to the T(4S) resonance at 
y/s = 10.6 GeV. However, due to the large width 
of the T(4S) (20 MeV), there is no significant enhance- 
ment to the continuum cross section from the resonance. 
Using the cut discussed above we obtain er eg _> 7 = 
(.g/lO^GcV -1 ) 2 1.2 10~ 5 pb, compared to a e ^ lv v = 
6.2 10~ 7 pb. 

Taking a total integrated luminosity for PEPII of 
312 /6 -1 , and cuts in pseudorapidity [—2, 2] and photon 
energy, the number of Standard Model events is almost 
zero. Requiring three <f> events gives a bound at 95% 
confidence: 

g < 8.9 10~ 6 GeV" 1 , (0.14) 

for m < Q.YlGeV. 

The integrated luminosity for KEKB at the time of 
writing is L = 5 10 5 pb . Imposing cuts as above in 
rapidity and in the photon energy the number of Stan- 
dard Model events is again almost zero. Requiring three 
events gives a bound at 95% confidence: 

g < 8.2 10~ 6 GeV" 1 , (0.15) 



for m < O.lSGeV. 

If luminosity increases to L — 1.0 10 6 pb -1 , almost 1 
Standard Model events is expected. Requiring 4 events 
gives a bound at 95% confidence: 

g < 5.9 10~ 6 GeV" 1 , (0.16) 

for m < 0A6GeV. 

For SupcrKEKB the expected luminosity per year is 
L = 1.0 10 7 pb -1 . In a two year period, imposing similar 
cuts, the expected number of background events is 12.4. 
Requiring 8.8 additional events gives a bound at 95% 
confidence: 

g < 1.9 10~ 6 GeV -1 , (0.17) 
for m<j, < 0.2SGeV. 
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